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ABSTRACT 

We study the large-scale clustering of galaxies in the overlap region of the Baryon 
Oscillation Spectroscopic Survey (BOSS) CMASS sample and the WiggleZ Dark 
Energy Survey. We calculate the auto-correlation and cross-correlation functions in 
the overlap region of the two datasets and detect a Baryon Acoustic Oscillation 
(BAO) signal in each of them. The BAO measurement from the cross-correlation 
function represents the first such detection between two different galaxy surveys. 
After applying density-field reconstruction we report distance-scale measurements 
Dyr^'^/rs = (1970±45, 2132±65, 2100±200) Mpc from CMASS, the cross-correlation 
and WiggleZ, respectively. The distance scales derived from the two datasets are con¬ 
sistent, and are also robust against switching the displacement fields used for recon¬ 
struction between the two surveys. We use correlated mock realizations to calculate 
the covariance between the three BAO constraints. This approach can be used to 
construct a correlation matrix, permitting for the first time a rigorous combination 
of WiggleZ and CMASS BAO measurements. Using a volume-scaling technique, our 
result can also be used to combine WiggleZ and future CMASS DR12 results. Finally, 
we show that the relative velocity effect, a possible source of systematic uncertainty 
for the BAO technique, is consistent with zero for our samples. 

Key -words: surveys, cosmology: observations, dark energy, distance scale, large scale 
structure of Universe 


1 INTRODUCTION 

The Baryon Acoustic Oscillation (BAO) signal is a relict of 


sound-waves became imprinted in the distribution of Cosmic 
Microwave Background (CMB) photons as well as in the 
matter density field. Over time, the density field evolved 
through gravitational collapse and cosmic expansion. While 
gravitational interaction can smear out the BAO signal, a 
complete destruction would require interactions over very 
large scales (today « 150Mpc), making the BAO feature a 
very robust observable. 


* E-mail: fbeutler@lbl.gov 


the early Universe, where photon pre ssure caused sound¬ 
waves to move out of over-densities ((Peebles fc Yul Il97d : 


ISunvaev fc Zeldovichll970l : lBond fc Efstathioulll987 l. These 


The BAO signal in the density field at different 
redshifts can be related to the BAO signal measured 
in the CMB and therefore allo ws employment of the so 
called standard ruler t echnique jBlake k, GlazebrookI 12003 : 
ISeo &: EisensteiDll2003h . We can compare the apparent size 
of the BAO signal measured in galaxy surveys with the 
absolute size of this signal measured in the CMB and use 
this to map out the expansion history of the Universe. 
Simulations have shown that the BAO signal is unaffected 


by systematic uncertainties down 

to the sub-percent 

level (lEisenstein. Seo & White 2007 

: Guzik & Bernstein! 

20071: Smith et al.l 2007, l2008l: 

Aneulo et al. !2008!: 

Padmanabhan & White! 20091: Mehta et al. 20 Ilf) and 


hence represents one of the most reliable tools available for 
precision cosmology. 


The most precise BAO measurement has recently been 
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reported by the Baryon Oscillation Spectroscop ic Survey 
(BOSS) collaboration (1%. lAnderson et al.ll2oi3 1 at a red- 
shift of z = 0.57. BOSS also achieved a 2% BAO distance 



reported a 4% BAO constraint. Future galaxy surveys point 
to Baryon Acoustic Oscillations as their main science driver 
(e.g. DESI. ISchekel et al.l[2oTll and will push well below 
percent-level precision over a large redshift range. 

In this paper we study the BAO signal in the overlap re¬ 
gion of the two largest galaxy sur veys available to date, the 
CMASS sample of BO SS DRll (lAh n et al. 1120141 ') and the 
WiggleZ galaxy survey jDrinkwater et ah 2010l 'l. We use the 
BOSS DRll sample, since the overlap between BOSS and 
WiggleZ does not increase with DR12. The two surveys have 
been produced independently, being subject to different see¬ 
ing conditions, instrumental noise, reduction pipelines and 
observation strategies. While BOSS targeted mainly Lumi¬ 
nous Red Galaxies (LRGs) and has a galaxy bias around 
6 « 2, WiggleZ selected blue Emission Line Galaxies (ELGs) 
with a bias around 6 « 1. Since we restrict our analysis to 
the overlap region between the two surveys, our datasets 
have the same sample variance. This allows a test of possible 
systematic uncertainties in galaxy clustering measurements. 
In this study we focus on BAO measurements while a com¬ 
panion paper ilMarin et al.ll2015l l investigates redshift-space 
distortions in the overlap region. 

The typical candidates for possible systematic errors in 
galaxy clustering studies are redshift space distortions as 
well as non-linear clustering, accompanied by a non-linear 
galaxy bias. Non-linear evolution is normally connected to 
the density peaks of the matter density field and hence we 
expect that galaxies which sit in high density regions, such 
as LRGs, should have stronger non-linear effects compared 
to field galaxies, such as ELGs. On the other hand, field 
galaxies have a larger redshift space distortion signal, which 
is a source of possible systematics if not modeled correctly. 
In addition to the properties of the density field itself, there 
are also survey specific aspects such as incompleteness that 
can modify the measured clustering statistic. 

Another possible systematic bias for BAO measure¬ 
ments is the relative velocity effect. The different velocities 
of baryons and dark matter after decoupling can affect early 
galaxy formation processes. At places where the relative ve¬ 
locity is large, baryons can escape gravitational potentials. 
The first stars in the Universe will form wherever the relative 
velocity is low, since this is where the baryons will condense 
into the gravitational potentials of the dark matter. The rel¬ 
ative velocity decays with redshift and is negligible at the 
redshift of the BOSS and WiggleZ galaxies. The question is 
whether the modulation of structure formation due to the 
relative velocity effect at high redshift is carrie d in the tracer 
galaxi es we observe today to measure BAO. IYoo fc SeliakI 
(l2013tl argue that galaxies which consist predominantly of 
old stars (such as LRGs) could carry this effect, while young 
galaxies might not. This would result in a modulation of the 
BAO signal in BOSS, but not in WiggleZ. In the overlap re¬ 


gion between BOSS and WiggleZ we can measure the BAO 
position in the same volume and compare whether the two 
surveys yield the same result, allowing the placement of con¬ 
straints on the possible impact of the relative velocity effect. 

In the next section we will introduce the two datasets 
used in our analysis, BOSS-CMASS and WiggleZ. Section[3] 
describes our correlation function estimate followed by a dis¬ 
cussion of our mock catalogues in Section [4] In Section [5] we 
present our technique of density field reconstruction followed 
by a discussion of our model for the correlation functions in 
Section!^ We then compare the obtained displacement fields 
and perform the correlation function fits in Section[7l In Sec¬ 
tion [H] we determine the correlation between BOSS-GMASS 
and WiggleZ. In Section[9]we introduce the relative velocity 
effect and perform fits to the data to constrain the relative 
velocity bias. We conclude in Section (TU] 

For clarity we will use the name GMASS-BW, WiggleZ- 
BW and cc-BW for the GMASS, WiggleZ and cross¬ 
correlation results limited to the overlap region between 
the two surveys. We adopt a flat ACDM cosmology with 
Dm = 0.27. The same model is used to construct templates 
for the BAO fits and hence our measurements should be 
used in conjunction with r^‘^{zd) ~ 150.18 MpcQ. 


2 DATASETS 


2.1 The BOSS survey 

BOSS, as part o f SDSS-III (lEisenstein et al.l I 2 OIII : 

I Dawson et al.l l2013l'l measured spectroscopic redshifts of 
« 1.5 million galaxies (and 150 000 q uasars) making use 
of the SDSS multi-fibre spectrographs llBolton et al.lf201^ : 
ISmee et al. I [ 2 ^. The galaxies are selected from multi¬ 
color SDSS imaging llFuku gita et al.ll 19961:iGunn et ^ Il998l : 
ISmith et al.l I 20 M : iGunn et al. 20061 : Doi et akf 201C 1 and 


cover a redshift range of z = 0.15 - 0.7, where the sur¬ 
vey is split into two samples called LOWZ {z = 0.15 - 0.43) 
and GMASS {z = 0.43 - 0.7). In this analysis we only use 
the GMASS sample. The survey is optimized for the mea¬ 
surement of the BAO scale and hence covers a large cos¬ 
mic volume (I4fr = 2.31 x 10®[Mpc/h,]®) with a density of 
n « 3 X 10“‘^[/i/Mpc]®, high enough to ensure that shot 
noise is not the d ominant error contribution at the BAO 
scale (IWhitell201llj . Most CMASS galaxies are red with a 
prominent 4000 A break in their spectral energy distribu¬ 
tion. Halo Occupation studies have shown that galaxies se¬ 
lected like the CMASS galaxies are mainly central galaxies 
residing in dark matter halos of IO^^Mq/Zi, with a 5 - 10% 
satellite fraction (IWhitell201lh . CMASS galaxies are highly 
biased (fo ~ 2), which boosts the clustering signal including 
BAO with respect to the shot noise level. 

We include three different incompleteness weights 
to account for sho rtco mings of the CMA SS dataset 
fsee lRoss et al.ir2012l and I Anderson et al.ll20l3 for details): 
A redshift failure weight, Wrf, a fibre collision weight, wtc, 
and a systematics weight, uigys, which is a combination of a 
stellar density weight and a seeing condition weight. Each 


1 


Sound horizon scale calculated with CAMB 


ll Lewis et al .l200(t) . 


© 0000 RAS, MNRAS 000, 000-000 









































































The BOSS-WiggleZ overlap region I 3 



RA n RA n 


Figure 1. Sky coverage of BOSS-CMASS DRll (black) and WiggleZ (red). The left plot shows the north galactic cap (NGC), while 
the right plot shows the south galactic cap (SGC). Five of the six WiggleZ regions are covered by CMASS, with region S22 being only 
partly covered. We only plot a random fraction of 3% of all galaxies. 


galaxy is thus counted as 

Wc = (Wrf + Wfc — l)Wsys. (1) 

Figure [T] shows the sky coverage of BOSS-CMASS with the 
north galactic cap (NGC) on the left and the south galactic 
cap (SGC) on the right. 

2.2 The WiggleZ survey 

The WiggleZ Dark Energy Survey dPrinkwater et al1l201Cll ~) 
is a large-scale galaxy redshift survey of bright emission¬ 
line galaxies (ELGs), which was carried out at the Anglo- 
Australian Telescope between August 2006 and Januar y 
2011 using the AAOmega spectrograph llSharp et al.ll200'^ . 
Targets were selected via joint ultraviolet and optical mag¬ 
nitude and color cuts using input imag ing from the Galax y 
Evolution Explorer (GALEX) satellite (iMartin et al.l[2005h . 
The survey is now complete, comprising 240 000 redshifts 
and covering 816 deg^ in six separate sky areas. The red- 
shift range is roughly 0.1 < « < 1.0 with a mean redshift at 
z = 0.6. Figure [T] shows the sky coverage of WiggleZ (red), 
where we name the different patches SI, S3, S22, N9, Nil, 
and N15. 

The strategy for completeness correction in WiggleZ is 
different to the method used in GMASS. Instead of weighting 
the data galaxies, the inco mpleteness has be en introduced 
into the random catalogues (iBlake et al.ll20l5 j and hence no 
completeness weighting is needed for this dataset. 

2.3 Definition of the overlap region 

We dehne the overlap region between CMASS and WiggleZ 
by splitting the sky into 0.1 deg^ bins and selecting all bins 
which contain CMASS as well as WiggleZ random galaxies. 
The redshift range is defined by CMASS and is given by 
0.43 < 2 < 0.7. Figure [T] shows the six WiggleZ regions 
(red), of which hve are covered by the BOSS-CMASS sample 
(black), with region S22 being only partly covered. 

The five overlap regions are shown separately in Fig¬ 
ure [2] Most of the incompleteness in these plots is caused 


Table 1. Effective volume and number of galaxies of the five 
distinct CMASS-WiggleZ overlap regions (see Figure as well 
as the total volume of the two surveys. The effective volume is 
calculated using Eq. [2] and Pq = 20 000/i“® Mpc® for CMASS 
and Pq = 5000/i“® Mpc® for WiggleZ. The names CMASS-BW 
and WiggleZ-BW stand for the CMASS and WiggleZ samples 
restricted to the overlap region between the two. 



Region 

VeS 

[lO'^h-^Mpc®] 

Agal 

SI 

CMASS-BW 

1.8 

5742 

WiggleZ-BW 

0.8 

6621 

S22 

CMASS-BW 

1.9 

6070 

WiggleZ-BW 

1.1 

10 339 

N9 

CMASS-BW 

3.1 

9356 

WiggleZ-BW 

1.7 

13 960 

Nil 

CMASS-BW 

3.7 

10 280 

WiggleZ-BW 

2.0 

15 324 

N15 

CMASS-BW 

4.6 

14 673 

WiggleZ-BW 

2.6 

22 736 

combined 

CMASS-BW 

15.1 

48 570 

WiggleZ-BW 

8.2 

71407 


CMASS-DR12 

232.2 

786 324 


CMASS-DRll 

204.0 

690827 


CMASS-DRIO 

150.6 

544 133 


CMASS-DR9 

76.9 

264 281 


WiggleZ 

21.3 

191 732 


by the WiggleZ survey, while the empty stripes in region 
Nil are a result of missing photometry and hence missing 
galaxies in CMASS. 

The relative importance of sample variance and shot 
noise in a galaxy clustering measurement is determined by 
the quantity n{z)P{k), where n{z) is the galaxy number den¬ 
sity and P{k) is the galaxy power spectrum amplitude at the 
BAO scale. Therefore we can trade a smaller galaxy density 
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Figure 2. The overlap region between BOSS-CMASS (black) and WiggleZ (red). Most of the angular incompleteness is caused by 
WiggleZ, while the empty stripes in region Nil are caused by incomplete photometric data in CMASS. To generate these regions, we 
divided the sky into 0.1 deg^ bins and included all bins which contain CMASS as well as WiggleZ random galaxies. We only plot a 
random fraction of 10% of all galaxies. 


with a larger galaxy bias and vice versa. The WiggleZ survey 
has a higher galaxy number density compared to CMASS, 
while CMASS galaxies have a larger bias and hence a larger 
power spectrum amplitude. The CMASS sample has been 
designed with the target nP = 3, while WiggleZ has aimed 
for nP = 1. The best quantity to compare the two surveys 
is the effective volu me, where we use the equation suggested 
by iTegmarkI (|l997h : 


Kff = 



n{x)Po 
1 + n{x)Po 


( 2 ) 


Here Pq is fixed to the amplitude of the power spectrum 
at the wavenumber of the first BAO peak, k « 0.06h/Mpc, 
which turns out to be Po = 20 000h“® Mpc^ for CMASS 
and Po = 5000h~^ Mpc^ for WiggleZ. The larger value of 
nP in the CMASS sample leads to a larger effective volume 
compared to WiggleZ (by about a factor of 2). The different 
volumes for CMASS and WiggleZ in each overlap region, as 
well as the combined volumes, are summarized in Table [T] 
The redshift distribution for the two samples limited to the 
overlap region is plotted in Figure El 


3 ESTIMATING THE CORRELATION 
EUNCTION 

We calculate the correlation function by counting the num¬ 
ber of galaxy-galaxy pairs, DD{s), as a function of scale 
s, as well as galaxy-random, DR{s), and random-random, 
RR{s) pairs. We then use the correlation function estimator 



Figure 3. Redshift distribution of CMASS-BW (red) and 
WiggleZ-BW (blue) combining the five separate regions. 


suggested by iLandv fc Szal^ (Il993l) : 


e(s) = i + 


DD{s) 

RR{s) 


- 2 


DR{s) 

RR{s) 

Na 


(3) 


where Ur = rid = ^i(®) represent 

the sums over the weights for all random and data galax- 
ie s, respectively. We include the inv erse density weighting 
of iFeldman. Kaiser fc PeacockI (1 19941 ): 


•Wi{x) = 


1 -I- n{x)Po ’ 


(4) 
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with Pq = 20 000/i“® Mpc® for CMASS and Po = 

5 000/i“® Mpc® for WiggleZ. In the case of CMASS we also 
include the completeness weighting of Eq. [T] 

The cross-correlation function between two tracers A 
and B can be calculated as 

^ D^D^{s) D^R^{s)(n^ 

- 1+ ^A^B(,) I^„a 

R^D^{s) hnf\ 

■ R^R^is) [T^J ■ 

Figure |4] shows the correlation functions for CMASS-BW 
(top), WiggleZ-BW (bottom) and the cross-correlation (cc- 
BW, middle). The three figures on the left present the results 
before density field reconstruction, while the figures on the 
right display the result after reconstruction (see Section [5] 
for a discussion of our reconstruction technique). The five 
correlation functions for the individual regions of CMASS- 
BW and WiggleZ-BW are indicated as grey lines. Using the 
covariance matrix (see next section) we can combine the 
correlation functions of the five sub-regions (colored data 
points). The auto-correlation functions of both surveys, as 
well as the cross-correlation function before and after recon¬ 
struction, show a clear BAO signal at around lOOMpc/h. 

We also measured the cross-correlation coefficient de- 
hned as 



^(s) = 


Ccc-BwC®) 


^CMASS-BW (s)^WiggleZ-BW (s) 


( 6 ) 


and presented in Figure [5] In linear theory we expect this 
quantity to be 


2 , , _ [^ + + h-a) + 

nheo..(^) - (1 + i/3A + i/3i)(l + i/lB + \Piy 

Assuming Ija = 2, fos = 1 and f = bP = 0.76 results in 
’'theory = 0.997. This expectation is included in Figure [5] 
(black dashed line). The mock realizations show a large cor¬ 
relation coefficient after density field reconstruction. We cur¬ 
rently do not have a model for the correlation function shape 
post reconstruction llWhitell2015l i and therefore we only use 
the pre-reconstruction result in our fitting in Section [5] 

Figure[5]also shows a small correlation coefficient before 
reconstruction (blue data points). We used Gaussian error 
distribution to obtain the uncertainties on the data points 
in Figure [5l however, the errors on r have a signihcant non- 
Gaussian component. In iMarin et al.l (l2015l l we performed 
fits to the auto- and cross-correlation functions having r as a 
free parameter. We find that r is consistent with 1 f or scales 
above 20Mpc//i (see Figure 5 in iMarin et aLllioisI l. 


4 MOCK REALISATIONS 


We produced 480 mock catalogues for each of 
the f ive overlap regions using the COLA tech¬ 
nique dTassev. Zaldarriaga fc Eisensteinll2013h . These mock 
catalogues will be presented in a separate paper together 
with t he details of the COLA implementation we em¬ 
ployed (iKoda et al. 2015|b These mock catalogues h a ve als o 
been used in Kazin et al.l ll20l4 l and iMarin et al.l ll2015h . 
Each simulation uses 1296® particles in a [600Mpc//i]® box 
resulting in a particle mass of 7.5 x IO^Mq/H, allowing 
us to resolve CMASS as well as WiggleZ size halos. The 






20 25 30 35 40 45 50 55 

s [Mpc/h] 


pre-recon data 
post-recon data 
pre-recon mock mean 
post-recon mock mean 


Figure 5. The measured correlation coefficient before (blue) and 
after (black) density field reconstruction. The dashed line shows 
the expectation of linear theory. The blue data points are shifted 
by 0.5 Mpc/h to the right for clarity. The solid lines indicate the 
mean correlation coefficient of the mock realizations before (red) 
and after (magenta) density field reconstruction. The error on 
the data points is derived from the variations in the 480 mock 
catalogues (grey lines). When fitting the correlation coefficient in 
Section|9]we only use the data before reconstruction, since we do 
not have a model for the correlation function post reconstruction. 


halos are identified using a friend-of-friend algorithm with a 
linking length of 0.2 times the mean particle separation. We 
use Halo Occupation distribution models to populate these 
halos with galaxies so that the mock realizations match 
the measured projected correlation functions Wp{rp), where 
Tp is the angular separation between a galaxy pair. The 
fiducial cosmology of these mock catalogues is flat ACDM 
with Qrn = 0.273, Iff, = 0.0456, Hq = 70.5 Mpc/km/s, 
(Ts = 0.812 and Us = 0.96. The comparison of the correlation 
functions measured in the mock catalogues and the data 
correlation functions are shown in Figure (6] The mocks 
match the WiggleZ and CMASS clustering on large scales 
while they slightly overestimate the clustering amplitude 
of the cross-correlation function. The discrepancies are less 
significant after reconstruction (three panels on the right in 
Figure [6]). 


4.1 Covariance matrix 

Using the mock realizations of the individual sub-regions 
we can produce covariance matrices for each of the auto- 
and cross-correlation functions. We calculate the covariance 
matrix using 

.. 480 

Cif = [^n(si) — C(s*)] “ S(®j)] ! (8) 

n=l 

with ^rt(si) being the n-th correlation function estimate at 
separation Si and the sum is over all 480 mock realizations. 
The mean value is defined as 

.. 480 

( 9 ) 
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Figure 4. The correlation functions of CMASS-BW (top), WiggleZ-BW (bottom) and the cross-correlation (middle) in the overlap 
regions between CMASS and WiggleZ. The grey lines show the correlation functions for the five individual sub-regions (see Figure [J, 
while the colored data points show the combined correlation functions calculated from Eq. 1101 The error bars are the diagonal of the 
combined covariance matrices (see Figure O. Note, that the scatter in the grey lines does not represent the error in the data points, 
since each grey line corresponds to a different volume and is weighted accordingly. The black lines show the best fit to the individual 
correlation functions corresponding to the upper part of Table [2] 


Instead of analyzing the 10 auto-correlation functions and 
five cross-correlation functions individually, we chose to 
combine the correlation functions to obtain two auto¬ 
correlation functions for CMASS-BW and WiggleZ-BW as 
well as one cross-correlation function. We combined the cor¬ 
relation functions of the five individual sub-regions using 
the covariance matric es calculated abo ve and follow i ng th e 
procedure outlined in I Whit j (l201ll) an d iBlake et al.l ll201ll) . 
Each sub-region is weighted by its corresponding uncertainty 

^ [cr'r(s), (10) 

regions i 


with C* being the covariance matrices of the individual sub- 
regions. The inverse covariance matrix for the combined cor¬ 
relation functions is given by 

c-i= ^ [cr\ (11) 

regions i 

which follows from Eq. 1101 The combined covariance ma¬ 
trices before and after density field reconstruction are pre¬ 
sented in Figure [3 The combined correlation functions for 
CMASS-BW, WiggleZ-BW and the cross-correlation func¬ 
tion are shown in Figure |4] as colored data points. We 
also compare the CMASS-BW correlation function with 
the CMASS-DRll correlation function in Figure [5] While 
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; cc-BW pre-reck>n 
■60 ^ 1 I I I I I I I I I I I 

0 20 40 60 80 100 120 140 160 180 200 

s [Mpc/h] 




0 20 40 60 80 100 120 140 160 180 200 

s [Mpc/h] 



Figure 6. Comparison between the mock realizations and the CMASS-BW (red, top) WiggleZ-BW (blue, bottom) and cross-correlations 
functions (black, middle). The grey lines show the 480 mock realizations with the mean given by the colored data points and the error 
representing the variance. The colored lines indicate the measurement in the data. The dashed black lines mark the fitting range which 
goes from 50 Mpc/h to 180 Mpc/h. 


the CMASS-BW correlation function before reconstruction 
is in excellent agreement with CMASS-DRll, we find the 
prominent BAO peak at slightly larger scales compared to 
CMASS-DRll. We will discuss this aspect further when fit¬ 
ting these correlation functions in Section [7.21 

It has been shown that the inverse covariance C~^ de¬ 
rived from a fi nite number of r e alizations underestimates the 
uncertainties llAndersonll200,ll : iHartlao. Simon fc Schneider! 
l2006l : [^rcival et al.ll2013l ). In the case of Gaussian errors and 
statistically independent bins, this effect can be accounted 


for by multiplying the variance estimated from the likelihood 
distribution with 


1 -b B{Nhins — Ap) 
""" - 1 -b 2A -b B{Np + 1) ’ 


( 12 ) 


where Amocks is the number of mock realizations, Abins is 
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8 Florian Beutler et al. 



Figure 7. The correlation matrix for the combined CMASS-BW, cross-correlation (cc-BW) and WiggleZ-BW correlation functions before 
(left) and after (right) density field reconstruction. These matrices are combinations of the individual matrices for the five separate regions 
using Eg. llll For each region we use 480 mock realizations for CMASS-BW and WiggleZ-BW. The color indicates the level of correlation, 
where red is high correlation, green is no correlation and blue is high anti-correlation. Since each set of CMASS-BW and WiggleZ-BW 
mock catalogues has been produced from the same simulation (see Section 0, there is a considerable amount of correlation between 
the three correlation functions, mimicking the situation of the real data. Given that we use a fitting range of r = 50 - 180Mpc/h with 
5Mpc//i bins, this matrix has 26 x 26 bins for each correlation function and 78 x 78 bins in total. 


ISOr 


- pre-recon 




50- 


0 - 


• CMASS-BW 

- CMASS-DR11 

I I I I I I I I I 


0 20 40 60 80 100 120 140 160 180 200 

s [Mpc/h] 



Figure 8. Comparison of the CMASS-BW (red data points) and CMASS-DRll (red dashed lines) correlation functions before (left) and 
after (right) density field reconstruction. Post reconstruction one can see that the BAO peak in CMASS-BW is at larger scales compared 
to CMASS-DRll, which leads to a smaller value of Dy as is also visible in the resulting likelihood distribution (see Figure fT^ . 


the number of bins, Np is the number of free parameters and 

1_ ^ 3 '! 

(-/Vmocks — A^bins — l)(A^mocks “ A^bins “ 4) ’ 

B = A(A^^ocks - A^bins - 2). (14) 


Furthermore, the sample variance needs to be multiplied by 


Nrr 


1 


A^mocks A^bins 2 


(15) 


Since the bins in a correlation function are not statistically 
independent, these correction factors are only an approxima¬ 
tion. Given the 480 mock realizations in our analysis, 26 bins 


and 5 free parameters (see Section[7(2}, we have ma = 1.033 
and = 1.095. However, when fitting all three correlation 
functions simultaneously (78 bins), this factor can rise to 
m„ = 1.4, significantly contributing to our error budget. 


5 DENSITY FIELD RECONSTRUCTION 

In linear theory, the co-moving position of the BAO peak is 
set after the epoch of decoupling, providing the foundation 
of its use as a standard ruler. There are, however, non-linear 
effects, which can change the BAO peak position, as well as 
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its shape. The mo st siRnificant effect is non-linear dampinR 
of the BAO peak (lEisenstein, Seo fc Whit3l2007l : ISeo et all 
l2008l i . This effect is often modeled with a Gaussian damping 
term. Although we are interested in the BAO peak position 
and not in its amplitude, damping can shift the peak posi¬ 
tion, because of the non-symmetri c shape of the correlation 
funct i on around the BAO peak jEisensteinjSeo&_WhiM 
20071: iGuzik fc Bernsteliil l2007l: ISmithet al l l2007l . l2008l: 


Amgulo_et_alJ_ 200 ^: iMefto e t ^ 1201 il l. Additionally, 
Crocce fc Scoccimarrol ( 2008l l found that mode coupling 


can lead to shifts in the BAO peak position. Interestingly, 
mode coupling as well as non-linear damping can be 
removed by a technique called density field reconstruction, 
meaning that the measured distribution of galaxies itself 
can be used to reduce the impact of these non-linear 
effects by estimating the displacements of galaxies from 
their i nitial position in the density field (|Eisenstein_eLalJ 
I2OO7I : iPadmanabhan fc White! l2009l : iMehta et al.l 2OIIII . 
Density field reconstruction enhances the signal-to-noise 
ratio of the BAO signature using extra information 
contained in the higher order cor relations of the galaxy 
distribution dEisenstein et al.1l2007l l. We apply density field 
reconstructio n to the observed d e nsity field following the 
formalism of iPadmanabhan et al.l (l2012li . First we smooth 
the observed and random fields with a Gaussian filter of 
the form 


G{k) = exp 


-(fcE, 


(16) 


where we choose Esmooth = 15Mpc//i (|Xu et al.ll2012l i. The 
over-density field is then calculated in real-space as 




(17) 


with Pg{x) and Prix) being the density of the smoothed 
galaxy and random distribution, respectively. The normal¬ 
ization is defined as 


nr _ Yjf’wijx) 

In linear perturbation theory, the displacement field 

T(a :) is related to _ the redshift-space density field 

by (jNusser fc DavislIlQQll i 

V-T(x)+/?V- [Tio,(f)] (19) 

where Tios is the line-of-sight component of the displace¬ 
ment field. Transforming this equation into Fourier space 
and using the approach ())(a;) = Efc ')’(^) exp(ifca;) and 5{x) = 
S{k) exp(ifea:) we get 

— (I>{k) [fc^ + ky + k^(l + /?)] = —(20) 

which we solve for <))(fc) for every wavenumber k. The 
displacement field is than given by 5'(fc) = —ij^(^(fc), 
which we Fourier transfer back into configuration space. 
Our approach uses the plane-parallel approximation, which 
is valid for the small angular coverage of the fiv e individual 
fields studied in this analysis llBlake et al.ll20llh . 


We then apply the displacement to our galaxies by shift¬ 


ing their x, y and line-of sight positions following 

Sios = Sold - (1 + /)^ios(:r), (21) 

Sa;,T/ — Sold ^ x,y{x) • (22) 


We do not apply the factor of (1 -f /) in the case of the ran¬ 
dom galaxies, since the random distribution does not contain 
redshift space distortions. During reconstruction we use the 
growth rate / = 0.7 as well as t he linear bias fo = 1.9 for 
GM ASS-BW (iBeutler et al.ll20l4l and 6 = 1.0 for WiggleZ- 
BW llBlake et al.l 2^^7'The three plots on the right of Fig¬ 
ure [J] show the correlation functions for GMASS-BW (top), 
WiggleZ-BW (bottom) and cross-correlation (cc-BW, mid¬ 
dle) after applying density field reconstruction. 


6 MODELING THE LARGE SCALE 
CORRELATION FUNCTION 


Our model for the galaxy correlation function follows the 
procedure oflAnderson et al.l ll2013li. The galaxv correlation 
function is given by 

^(s) = P^^m(as) - 1 - A(s) 

(23) 

where 


,1 / \ 0,1 <22 

^(s) = “T H-h 03 . 

5^ S 

(24) 

The matter correlation function is 

through (lEisenstein. Seo & Whitell2007li 

obtained 

Cm(s) - J 2^2 P{k)jo{ks)e 

(25) 


with (Ts = 2 Mpc//i and the monopole power spectrum is 
given by 

P{k) = Psm,lin(fc) [l -f (Olin(fc) - 1) e" /"] . (26) 

We fix Eni = 8 . 8 Mpc/h, before reconstruction and 

Eni = 4.4 Mpc/fe aft er reconstruction ijAnderson et al.ll201.‘ll : 
iMagana et al.l I^Olll f. The function Oiin(fc) represents the 
oscillatory part of the fiducial linear power spectrum and 
Psm(fc) is the smooth power spectrum monopole. To obtain 
Psm.iin( fc) we fit the fiducial lin ear power spectrum, Piin(fc), 
with an lEisenstein fc Hul (Il998l l no-Wiggle power spectrum, 
Pnw(fc), together with five polynomial terms: 

Psm,lin(fc) = B Pnw(fc) -f ^ -|- C4 -|- Csfc. (27) 

The oscillatory part of the power spectrum is given by 


Olin(fc) 


Plin(fc) 

Psm,lin(^) 


(28) 


Our model in Eq. m has five free parameters (P, a, ai-s). 

To turn the constraint on a into a physical parameter 
we use 


DtrH^)rs 


(29) 


with 


Dv{z) 


(1 - 1 - z)^Da{z) 


cz 


1/3 


H{z) 


(30) 
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Figure 9. Comparison of the displacements in the CMASS-BW 
catalogue using the CMASS-BW displacement field (x-axis) and 
the WiggleZ-BW displacement field (y-axis). The plot shows a 
random selection of 5% of all galaxies in the five overlap regions. 


where Da{z) is the angular diameter distance and H{z) is 
the Hubble parameter. 


7 TESTING FOR BAO SYSTEMATICS 

Although the linear bias model was always believed to 
be sufficient for scales as large as the BAO signal, some 
studies using halo catalogues from N-body simulations 
suggest that the re are scale-dep e ndent bias effects even 
on BAO scales (iNoh et ahl l2009l : iDesiacaues et al.l l2010l : 
IWang fc Zhanll2013h . This means that the BAO signal can 
vary, depending on the tracer chosen to map the underly¬ 
ing density field. In the following sections of this paper, we 
will fit the large scale correlation function of CMASS and 
WiggleZ, and compare the displacement fields derived from 
the two surveys. Since we restrict our analysis to a common 
volume, we expect the results to be correlated. Since the two 
surveys trace the underlying density field differently, we can 
test for systematic effects in the BAO analysis. 


7.1 Comparing the CMASS and WiggleZ 
displacement fields 

In Section [5] we derived two displacement fields using the 
CMASS-BW and WiggleZ-BW galaxies, respectively. Here 
we are interested to learn (1) whether one of the displace¬ 
ment fields leads to better BAO constraints and (2) whether 
there are any systematic shifts in the BAO position depend¬ 
ing on which displacement field is used for the reconstruc¬ 
tion. 

We apply the displacement field derived using the Wig¬ 
gleZ survey to the CMASS galaxies and the displacement 
field derived from the CMASS survey to the WiggleZ galax¬ 
ies resulting in four datasets: 


1 CMASS-BW using the CMASS-BW displacement field, 

2 CMASS-BW using the WiggleZ-BW displacement field, 

3 WiggleZ-BW using the CMASS-BW displacement field, 

4 WiggleZ-BW using the WiggleZ-BW displacement field. 


Figure [9] compares the two displacements for each CMASS- 
BW galaxy. We quantify the correlation between the dis¬ 
placement fields using the correlation coefficient 


r(A,B) 





(31) 


where represents the mean of the displacement and = 
^^{xi — T)^, summing over all galaxies i. The correlation 
coefficient for the regions (SI, S22, N9, Nil, N15) = (0.65, 
0.64, 0.67, 0.68, 0.81) for the CMASS-BW galaxies and 
(0.49, 0.74, 0.76, 0.73, 0.75) for the WiggleZ-BW galaxies. 
The smallest region (SI) shows the lowest correlation coef¬ 
ficient, indicating that volume effects do play a role in this 
case. 

The mean difference between the CMASS-BW and 
WiggleZ-BW displacement fields using CMASS galaxies 

is = ^CMASS-BW _ ^WiggleZ-BW _ _q ^ 

0.016 Mpc//i, while for the WiggleZ catalogue = 

—0.150± 0.020 Mpc/h (the errors are the error on the mean 
between all galaxies). We therefore find moderate differ¬ 
ences between the two displacement fields. The difference 
between the two displacement fields does depend linearly on 
the amplitude of the displacement, (^^cmass + 4/wiggiez)/2, 
as shown in FigurefTOlfor the CMASS-BW galaxies (red) and 
the WiggleZ-BW galaxies (blue). Such a discrepancy could 
be caused by an incorrect assumption of the bias parameter 
when deriving the displacement field. 

To further investigate the impact of the two different 
displacement field on the BAO scale, we now calculate the 
correlation functions using both displacement fields. The 
correlation functions after combining the five different re¬ 
gions are presented in Figure [TT] In the next section we will 
fit these correlation functions and derive BAO constraints. 


7.2 Fitting the large scale correlation function 

We start with fitting the individual correlation functions 
of CMASS-BW and WiggleZ-BW as well as the cross¬ 
correlation function. We search for the best fitting parame¬ 
ters defined by the minimum given by 

= ( 32 ) 

ij 

where D is a vector containing the difference between 
the data and the model. Using the fitting range 50 — 
180 Mpc/h in 5 Mpc/h bins results in 26 elements for the 
vect or D. We use the python-ba sed MCMC sampler em¬ 
cee llForeman-Mackev et al.ll2013f l to derive the likelihood. 
The results are shown in Table [2] We can clearly see that 
the constraints for all three correlation functions improve 
significantly after reconstruction. The resulting BAO con¬ 
straints are worse, however, if we switch the displacement 
fields between the two surveys (these results are labeled as 
“(switched)” in Table[2)). The same result occurs in the mock 
realizations, where 58% of the WiggleZ-BW mock catalogues 
show a larger uncertainty on a when using the CMASS-BW 
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Figure 10. Comparison of the CMASS-BW (red) and WiggleZ- 
BW (blue) galaxies using the CMASS-BW and WiggleZ-BW dis¬ 
placement fields, respectively. The a;-axis is the mean displace¬ 
ment, while the j/-axis shows the difference. The CMASS-BW 
data points are shifted by 0.25 Mpc/h to the right for clarity. The 
uncertainties are derived from the mock realizations. 


displacement field for reconstruction instead of the WiggleZ- 
BW displacement field. Similarly 59% of the CMASS-BW 
mock realizations show poorer constraints when using the 
WiggleZ-BW displacement field. The resulting likelihood 
distributions for all fits are presented in Figure [12] including 
a comparison to the CMASS-DRll result. The likelihood 
distributions are reasonably approximated by Gaussians. 

In the limit of sample variance dominated uncertainties 
the three correlation functions would carry the same amount 
of information and only one of them would need to be ana¬ 
lyzed. In the shot noise limit all three correlation functions 
would be independent and would need to be analyzed to¬ 
gether to make maximal use of the available information. In 
the case of CMASS-BW and WiggleZ-BW, shot noise does 
contribute significantly to the error budget, so that a com¬ 
bined analysis is beneficial. 

Therefore we now fit all three correlation functions to¬ 
gether using the combined covariance matrix shown in Fig¬ 
ure (3 The fit has 13 free parameters, one scaling parame¬ 
ter, a, the three bias parameters, Bcmass-bw, Bcc-bw and 
fJwiggieZ-BW as well as three polynomial terms per correla¬ 
tion function. The data vector for this fit is given by 


^CMAp-BW(^^) _ ^CMASS-BW(^^) . 
tCMASS-BWz„ \ cCMASS-BW/„ \ 

Smodel l^26j “ ?data 


D = 


Ccc-Bww / ^ \ tcc-BW/'_ \ 

Smodel l^26j — Sdata (526j- 


^WiggleZ-BW 
Smodel 


(si)-Cif 


WiggleZ-BW 


(si) 


V 


^WiggleZ-BW/ \ _ ^WiggleZ-BW 
Smodel V Sdata 


(S26) 


/ 


(33) 





Figure 11. Comparison of the correlation functions for CMASS- 
BW (top), WiggleZ-BW (bottom) and the cross-correlation func¬ 
tion (cc-BW, middle) after density field reconstruction using 
the two displacement fields derived from the CMASS-BW and 
WiggleZ-BW galaxies, respectively. The label “(switched)” indi¬ 
cates the case where the displacement field of the other survey 
has been used for reconstruction. The red data points are shifted 
to the right by IMpc/h for clarity. 
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12 Florian Beutler et al. 


containing 3 x 26 = 78 bins in total. The results are 
shown in the lower part of Table [2] In the case of pre¬ 
reconstruction, this fit is driven by the cross-correlation 
function, which has significantly smaller uncertainties than 
any of the auto-correlation functions and leads to a value of 
a = 1.095 ± 0.068. After reconstruction it is the CMASS- 
BW constraint which drives the combined fit, leading to 
a = 0.966 ± 0.031. The combined constraint on a is worse 
than the CMASS-BW only constraint on a, which is mainly 
caused by the large scaling factor of = 1.4 (see Eq. [IS), 
needed for this fit. 


7.3 Comparison to mock realizations 


The question now is whether the measured a values for 
CMASS-BW and WiggleZ-BW are consistent. We can test 
this, by using the 480 correlated mock realizations, which we 
used to calculate the covariance matrix in Section SU We 
calculate the correlation function for each mock catalogue 
and repeat the fitting procedure described in the last sec¬ 
tion. Figure [TT] shows the distribution of the different con¬ 
straints on a for the mock realizations, together with the 
results found for the data (red data points). The ellipses in 
these plots are the Icr standard deviation including the cor¬ 
relation between the different measurements. We only plot 
results which have a value of a between 0.6 < a < 1.4 as 
well as an error on a less than 25%. For CMASS-BW we 
have 318 out of 480 mock catalogues which fulfill these cri¬ 
teria, while for WiggleZ-BW there are 242 and 302 for the 
cross-correlation. 

Using the mock realizations the standard deviations for 
Q are = (0.039, 0.036, 0.055) for CMASS-BW, cc-BW and 
WiggleZ-BW, respectively. While the mock realizations pre¬ 
dict the best BAO constraint to be in the cross-correlation 
function, in the data the most accurate distance scale mea¬ 
surement (post-reconstruction) is in CMASS-BW. This re¬ 
sult is, however, consistent with sample variance, and we 
have many mock realizations which show a similar behav¬ 
ior. Figure [TJ] compares the distribution of errors for the 
mock realizations with the data. The signal-to-noise ratio of 
the auto-correlation functions is given by 


uaPaa 
^aPaa + 1 


(34) 


while the signal- to-noise ratio of th e cross-correlation func¬ 
tion scales with (ISmith et al.ll200^ 1 


o ynAWnPAB 

\J{nAPAA + l){nBPBB + 1 ) + nAnsP^B 

Assuming nP = 3 for CMASS-BW and nP = 1 for WiggleZ- 
BW as well as Pcmass-bw = 4PwiggieZ-BW (Pcc-bw = 
2Pwiggiez-Bw), we find 0.75, 1.04 and 0.5 for CMASS-BW, 
cc-BW and WiggleZ-BW, respectively. This rough estimate 
of the expected signal-to-noise ratio agrees well with the 
result of the mock realizations. We can therefore conclude 
that in our actual data measurements the WiggleZ-BW con¬ 
straint has an error slightly larger than expected, the cross¬ 
correlation function error is close to the expected error and 
the CMASS-BW error is smaller than expected. 

To determine whether the three different BAO con¬ 
straints agree, we construct the vector T = (ctcmass-bw — 
Q!cc-bw, Qcmass-bw — awiggieZ-Bw)- Using the mock cat- 





CMASS-BW (red), WiggleZ-BW (blue) and the cross-correlation 
function (black), before (top) and after (middle) density field re- 
constructi on. The dashed red lin e shows the CMASS-DRll con¬ 
straint of llAnderson et al]|2013l l. The bottom panel displays the 
result where the displacement fields derived from the two surveys 
have been switched. 
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The BOSS-WiggleZ overlap region I 13 


Table 2. Summary of the fitting results. The first sector reports the fits to CMASS-BW, WiggleZ-BW and cross-correlation (cc-BW) 
functions individually, while the second sector shows the combined fits to all three correlation functions. For each case we list the result 
pre- and post reconstruction, as well as the result where we switched the displacement fields used for reconstruction (switched). 

The errors on each parameter are obtained by marginalizing over all other parameters. The fitting range is 50 - 180Mpc//i in 5Mpc//i 
bins leading to 26 bins and 5 free parameter (bias, three polynomials and a) in case of the fit to the individual correlation functions. 
When fitting all correlation functions simultaneously (last three rows) there are 3 X 26 = 78 degrees of freedom and 13 free parameters. 
The likelihood distributions are shown in FigureOur fiducial sound horizon is r®‘^( 2 ^) = 150.18 Mpc. 


survey 

OL 

[Mpc] 



CMASS-BW pre-recon 

i.uzy_o.o85 

2100Zlro 

31.2/(26 

-5) 

CMASS-BW post-recon 

0.970 ±0.022 

1970 ± 45 

22.6/(26 

-5) 

CMASS-BW post-recon (switched) 

0.976 ± 0.029 

1982 ± 59 

23.1/(26 

-5) 

cc-BW pre-recon 

1 nyo+O.Ubb 
l.UfO_Q Qg7 

21801““ 

22.8/(26 

-5) 

cc-BW post-recon 

1.050 ±0.032 

2132 ± 65 

26.7/(26 

-5) 

cc-BW post-recon (switched) 

1 noo+0.u5y 
J..UZO_Q Qg2 

2080l“0 

13.0/(26 

-5) 

WiggleZ-BW pre-recon 

I'OSto.ls 

219012« 

18.1/(26 

-5) 

WiggleZ-BW post-recon 

1.033 ±0.10 

2100 ± 200 

10.4/(26 

-5) 

WiggleZ-BW post-recon (switched) 


2190lg“ 

15.6/(26 

-5) 

combined-BW pre-recon 

1.095 ± 0.068 

2220 ± 140 

84.6/(78 - 

-13) 

combined-BW post-recon 

0.966 ± 0.031 

1956 ± 63 

103.5/(78 

-13) 

combined-BW post-recon (switched) 

u.y ^Z_Q 078 

1974l?|8 

60.5/(78 - 

-13) 



Figure 13. The error on a measured in the 480 mock realiza¬ 
tions of CMASS-BW (red), WiggleZ-BW (blue) and the cross¬ 
correlation function (black). The uncertainty CTc, is taken to be the 
mean of the upper and lower 68% confidence levels. The values 
measured in the data are shown by the correspondingly colored 
vertical dashed lines. 


alogues we obtain a correlation coefficient for the two com¬ 
ponents of the vector T of 60.6%, which allows us to con¬ 
struct a covariance matrix and calculate ~ T"'"C~^T 
with two degrees of freedom. Post reconstruction we get 
X^/d.o.f. = 1.25/2 when using the mean uncertainty pre¬ 
dicted from the mocks, and 4.1/2 when using (averaged) la 
errors from the actual measurements. The three different 
BAO measurements are therefore statistically consistent. 

The BAO peak in CMASS-BW leads to a distance con¬ 
straint of Dvr^'^/rg = 1970 ± 45 Mpc (post reconstruc¬ 
tion), which is low compared to the CMASS-DRll result of 


Dvrg‘^/ra = 2068 ± 20Mpc llAnderson et al.l [201311^ . Since 
CMASS-BW uses about 7% of the CMASS-DRll galax¬ 
ies, we can approximate the correlation coefficient, r^, to 
be about 0.07, implying that these constraints are fairly in¬ 
dependent. We can quantify the significance of the deviation 
between two correlated measurements with 


/—2-;-2-n-7-7 

\Jo\ + ag — 2x r((JA, aBjaAaB 

where r{aA,aB) is the correlation coefficient. This leads to a 
significance of 2.1cr for the deviation between CMASS-DRll 
and CMASS-BW. 

The WiggleZ survey has a distance constraint of 
Dvr^'^/rg = 2100 ± 200 Mpc, which we can compare t o 
Dvr^'^/vg = 2221 ^®q4 Mpc measured in iKazin et al.l ll2014h . 
The WiggleZ constraint has a slightly different redshift 
range compared to our WiggleZ-BW constraint (0.4 < 2 < 
0.8) and does include the additional sky region S3 (see Fig¬ 
ure [T]), therefore a direct comparison is not possible. 

IRoss et ll2014h split the CMASS sample based on 
fc-be corrected i-band absolute magnitudes and [r—ijo.ss col¬ 
ors, yielding two sub-samples with bias b = 1.65 (blue) and 
b = 2.3 (red). Studies of the BAO scale in these sub-samples 
revealed no statistically signifi cant deviati o ns, in agreement 
with our findings. Different to IRoss et ffi] (|2014| ) our study 
is based on two entirely different surveys and therefore also 
includes possible systematics due to instrumentation, tele¬ 
scope site conditions or reduction pipeline. 


^ This value has been corrected by the ratio of the fiducial sound 
horizons to allow direct comparison with our constraint. 
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Figure 14. These plots show the distribution of a for the mock 
realizations of CMASS-BW, WiggleZ-BW and cc-BW. We only 
plot results which have a value of a between 0.6 < a < 1.4 as 
well as an error on o less than 25%. The ellipse represents the 
1(T distribution drawn from the variance and correlation coeffi¬ 
cient. Note that the ellipse has not been derived from the black 
points drawn in these plots but instead from jack-knife samples to 
avoid outliers. The red datapoint shows our measurement post¬ 
reconstruction from Table[2] The agreement between the red data 
points and the black ellipse is only required if the cosmology of 
the mocks is the true cosmology. 


8 DETERMINING THE CORRELATION 

BETWEEN BAO CONSTRAINTS IN CMASS 
AND WIGGLEZ 

j.fid 

In the last section we set constraints on the quantity Dv -f— 
using the CMASS-BW and WiggleZ-BW auto-correlation 
functions as well as their cross-correlation function. The con¬ 
straint from the WiggleZ-BW auto-correlation function uses 
almost all WiggleZ information within the redshift range 
0.43 < z < 0.7, excluding only region S3, which is small in 
comparison. The entire CMASS sample, however, covers a 
sky area much larger than the overlap region. In this sec¬ 
tion we will determine the correlation of WiggleZ-BW and 
cross-correlation constraints fo und in the last s e ction with 
the CMASS-DRll constraint of lAnderson et ahl (l2013l ). We 
will then construct a covariance matrix which allows us to 
use our results together with the result of CMASS-DRll 
for cosmological parameter constraints. We will also provide 
an estimate of the correlation for the future CMASS-DR12 
constraint. 

First we divide the two surveys into the following sub- 
regions: 


1 BOSS-CMASS, excluding overlap region. 

2 BOSS-CMASS, in overlap region (CMASS-BW). 

3 WiggleZ, excluding overlap region. 

4 WiggleZ, in overlap region] (WiggleZ-BW). 


For each region, the parameter of interest is the constraint on 
the scaling parameter a. Assuming no correlation between 
the different regions, the final constraints for each survey are 
given by 


qb = 1 + 
aw = 1 + 


Vi(a?i — 1) -b V2 (q:2 — 1) 
V 1 + V 2 

Vs^as — 1 ) + 14 ( 0:4 — 1 ) 

14 + 14 


(37) 

(38) 


where the subscripted numbers refer to the four survey sub- 
regions described above, and V is the volume given in Ta¬ 
ble [T] The subscripted B stands for BOSS-CMASS and W 
stands for WiggleZ. The correlation between ob and ow is 
given by 


r (qb, Ow) 


r^(Q 2 ,04)1414 

1b Iw ’ 


(39) 


where the correlation coefficient between the constraints in 
the overlap regions, r^(o 2 , 04 ), can be calculated from the 
480 mock realizations as 


r^(02,04)= , E^^°(oz-02)(04-04) ^ 

7E"®°(«2-02)2E"®°(«4-04) = 


For practical reasons we use a jack-knife approach, in which 
we determine o for the mean of A — 1 realizations, excluding 
each of the 480 mock realizations in turn. The correlation 
coefficients between CMASS-BW and WiggleZ-BW in the 
overlap region is r^(o 2 , 04 ) = 0.301. Similarly, we can define 
the correlation coefficient between the auto-correlation func¬ 
tions and the cross-correlation function giving F{ac, 04 ) = 
0.570 and F{a 2 ,ac) = 0.584, where the subscripted 
C stands for the cross-correlation function. To obtain 
the correlation coefficient between our W iggleZ-BW con¬ 
strain t and the CMASS DRll constraint of lAnderson et al] 
1 I 2 OI 3 II . we use the volume of CMASS-DRll (as given in 
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Table [T]) in Eq. [39] and set V 4 = Vw- Therefore the 
data vector D = (CMASS-DRll, cc-BW, WiggleZ-BW) = 
(2056, 2132, 2100) Mpc has the following correlation matrix: 

/ 1 

Rdrii = I r^{aB, 

\ 

( 1 
= 0.043 

\ 0.022 

The covariance matrix is 
the vector V contains the variance of the individual con¬ 
straints. In our case we have Vdrii = (20,65,200) Mpc, 
wher e we adopted t he C MASS-DRll uncertainty (left) 
from [Anderson et al.l (I2OI3I ') together with the WiggleZ-BW 
(right) and cross-correlation function (middle) uncertainties 
from Table (21 This approach leads to the following covari¬ 
ance matrix 

/ 400 
Cdrii = j 56 
\ 88 

and the inverse is given by 
/ 250.47 

Cdrii = -3-48 35.11 1 x 10"®. (43) 

\ 0.09 -6.50 3.70 ) 


4225 (42) 

7410 40000 / 


ac) 1 

. 04 ) r^{ac,a4) 1 


0.570 1 


(41) 


given by C = V^RV, where 


Since the overlap volume between CMASS and WiggleZ 
will not change with DR12, this formalism can be rescaled 
to obtain the correlation between our results and future 
CMASS data releases. For example, using a cosmic volume 
of 2.322 /i“®Gpc® for DR12 we find the following correlation 
matrix 

Rdri2 = ( 0.038 1 ) . (44) 

V 0.020 0.570 1 / 


This covariance matrix is only correct assuming that any 
correlation between these surveys can be scaled with volume. 


So far we have only used WiggleZ galaxies in the redshift 
range 0.43 < z < 0.7, ignoring a significant fraction of Wig¬ 
gleZ galaxies at higher and lower redshifts. We can combine 
our results with the high redshift (0.6 < z < 1.0) WiggleZ 
measurement reported in iKazin et al.l (l2014h and given by 
Dvt^’^/ r 3 = 2516 ±86 Mpc. This measurement has an effec¬ 
tive redshift of Zeft = 0.73. The effective volume of WiggleZ 
in the redshift range 0.6 < z < 0.7 (overlap between the 
high redshift WiggleZ measurement and the CMASS red¬ 
shift range) is 6.1 x 10^/i“®Mpc®. The correlation matrix 
including the high redshift WiggleZ data point (labeled by 
ext.) would be 


Rr 


( 1 
0.043 
0.022 
\ 0.017 


1 

0.57 

0.42 


with the data vector 


y~^ext _ 


1 

0.51 


(45) 


1 / 


/CMASS-DR11\ 


/2056\ 

cc-BW 


2132 

WiggleZ-BW 


2100 

\ WiggleZ-highz y 


1^2516/ 


Mpc (46) 


and the variance vector is Rdrh = (20, 65, 200, 86) Mpc. 
This result makes the additional assumption that the cor¬ 
relation coefficient we found for the CMASS redshift range 
can be scaled to the overlap redshift range 0.6 < z < 0.7. 
We do not combine our measure ments with t he lo w and 
medium redshift bins reported in iKazin et al.l (l2014h since 
both overlap with the BOSS-LOWZ redshift range. There¬ 
fore the results reporte d in this section can be combined 
with the BOSS-LOWZ ([Anderson et al.ll2013l : lToieirdl2014[l 
measurement straightforwardly. 

The assumed hducial cosmologies used in the different 
measurements above are not the same, resulting in differ¬ 
ent fiducial sound horizons. The sound horizon used in the 
CMASS-DRll analysis is = 149.28 Mpc. IKazin et al.l 
|[2014[) have rf'* = 148.6 Mpc and our analysis uses = 
150.18 Mpc. When comparing the measurements of Dy 
above with a cosmological model one has to include the ratio 
of the hducial sound horizon and the sound horizon of that 
model. 


9 THE RELATIVE VELOCITY EFFECT 

In this section we discuss one possible source of systematic 
unce rtainty for BAO constraints, the relative velocity ef¬ 
fect (iTseliakhovich fc Hirata[l2010[ '). We hrst introduce the 
idea of the relative velocity effect and discuss our model, 
before htting the model to the data. 

While dark matter perturbations start to grow di¬ 
rectly after the end of inhation, baryon perturbations cannot 
grow until they decouple from the photons, about 380 000 
years later. The different velocities of dark matter and 
baryons after decoupling means th at there is a relative ve¬ 
locity between the two com ponents (ITseliakhovich fc Hiratal 
[201(]l : iFialkov et al.1 [20141 1 . The relative velocity can allow 
baryons to escape the dark matter potentials and prevent 
the formatio n of the hrst stars in regions with high rela¬ 
tive velocity dFi alkov et al.||2012|: [McQuinn fc O’LearvlIioi^ : 
iNaoz. Yoshida fc Gnedin[|2012[ ~). This modulation would se¬ 
lect regions with small relative velocity to hrst undergo 
reionisation. Since the relative velocity ehect decays with 
1 ± z, it mainly affects the high-redshift Universe. However, 
it has been speculated that galaxies which form at high 
redshift carry this selection process down to low redshift, 
perhaps through processes such as altering the metal abun - 
dances or supernovae feedback (IYoo. Dalai fc Seliak[|201l[ L 
In Fourier space the relative velocity ehect has an oscillatory 
pattern on large scales which is out of phas e with the Baryon 
Acoustic Oscillations (IYoo fc Seliakll^OlSlL In conhguration 
space the relative velocity ehect modihes the clustering am¬ 
plitude primarily below the sound horizon, lea ding to a shift 
of the BAO peak (ISleoian fc Eisensteinllioi^ L 

The hypothesis is that old galaxies still carry the selec¬ 
tion of the relative velocity ehect, while young galaxies do 
not. Under this hypothesis we can measure the relative ve¬ 
locity bias by comparing clustering statistics of BOSS and 
WiggleZ, since BOSS mainly targeted (old) luminous red 
galaxies, while WiggleZ selected young star forming galaxies 
(ELGs). Our analysis method therefore assumes that BOSS 
galaxies are ahected by the relative velocity ehect, while 
WiggleZ galaxies are not. 
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Figure 15. Comparison of the auto-correlation functions for different values of the velocity bias parameter. The plot on the left shows 
the entire correlation function, while the plot on the right focuses on the BAO peak. The relative velocity effect causes an increase in the 
correlation function amplitude as well as a shift in the BAO peak position towards smaller scales. For these figures we assumed bi = 2 
and b 2 = —0.4. The dashed lines represent the correction terms for different values of the velocity bias; the dashed lines added to the 
black solid line results in the correspondingly colored solid lines. 




Figure 16. Comparison of the cross-correlation functions for different values of the velocity bias parameter. The plot on the left shows 
the entire correlation function, while the plot on the right focuses on the BAO peak. The relative velocity effect causes a small shift 
in the BAO peak position towards smaller scales. Unlike for the auto-correlation function, the relative velocity effect does not change 
the amplitude of the cross-correlation function. For these plots we assumed 6^ = 2, 6^ = —0.4, bf = 1 and b^ = 1. The dashed 
lines represent the correction terms for different values of the velocity bias; the dashed lines added to the black solid line results in the 
correspondingly colored solid lines. 


9.1 Modeling 


To model the correlation function including t he relative ve¬ 
locity effect we follow the implementation of lYoo &; SeliakI 
(|2Q13l) . In this model, the galaxy density field is given by 

5g{x) = bi5m{x) + y [sl,{x) - a^] + b„ [Ur - CTur] , (47) 

where the relative velocity Ur is computed at the linear order 
and the matter density is computed to the second order. 
The auto-power spectrum from such a density field can be 


written as 

Pg(fc) = fe?PNL(fc) + I ■0^Pm{q)Pm{\k-q\) 

-I- 26ifo2P2(g, k - q) + 

+ 4&i6„P2(g, k - q)Gu{q, k - q) 

+ 2 b 2 bvGuiq, k- q) + 2blGl{q, k - g)j 
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with the kernels 


Gu(fci, fe) 


F 2 {k\, te) 


Tur(kl) Tur(k2) fcl • 


Tm{kl) Tm{k 2 ) klk 2 
2 

ki ■ k 2 \ ^ ^1 ■ 


kik2 


k 2 


(49) 



The dimensionless relative velocity transfer function Tu^ is 
defined as 


= 


Tv. — Tv. 


O' Vj- 


(51) 


where and Tv^^^ are the velocity transfer functions of 
baryons and cold dark matter, respectively and the normal¬ 
ization is given by 


, , 1 7 dfc n 


k 


0.002 Mpc 


(52) 


The transfer functions describe the evolution of each mode 
with redshift as 


S{k, z) = T{k, 2)5pri(fc), (53) 


where 5pri is the primordial density perturbation. Before 
reionisation the linear relative velocity transfer function does 
not change its shape but only its amplitude (see Figure fTTl) . 
Assuming the relative velocity effect has been imprinted be¬ 
fore reionisation, any redshift above reionisation can be cho¬ 
sen to calculate the relative velocity transfer functior0. We 
calculate the velocity transfer function T^ at z = 15, while 
the matter transfer function in the denominator of the ve¬ 
locity kernel is calculated at 2 = 0.57. 

We can turn this model into configuration space using 
a Fourier transform as given in Eq. 1251 Figure [TJ] shows 
the correlation function for different values of the relative 
velocity bias by . The relative velocity bias causes an increase 
in the correlation function amplitude as well as a shift of the 
BAO peak position towards smaller scales. 

Using the definition of Eq.jd^we can write a model for 
the cross-correlation function of a tracer A, which carries 
the relative velocity effect, with another tracer B, which 
does not carry this effect: 

P^^ik) = btbfP^^k) + I ^P^(g)P^(|fc - g|) 

+{bib2 + b2bf)F2{q,k — q)+ 

+ 2 bfb^F 2 {q,k-q)Gu{q,k-q) 

+ b2bvGu{q,k-q)]- 

Figure [TS] shows the cross-correlation function for differ¬ 
ent values of the relative velocity bias by As was the case 
for the auto-correlation function, the relative velocity bias 
causes a shift of the BAO peak position towards smaller 
scales. In the cross-correlation function the shift is about 
1/3 of the shift present in the auto-correlation function. Un¬ 
like the auto-correlation, the relative velocity effect does not 


® CMBfast and early versions of CAMS do not include the effect 
of reionisation to th e baryon transfer function (see Section 2.3 
in lLesgourguesll201 jf . 



Figure 17. The velocity transfer function Tv — Tv^.^^ rela¬ 

tive to the cold dark matter transfer function Tcdm ■ While at high 
redshift this ratio only changes in amplitude, at low redshift it has 
additional contributions from reionisation. Assuming that the rel¬ 
ative velocity effect has been imprinted before reionisation, the 
velocity transfer function can be evaluated at any redshift (above 
reionisation), since the amplitude is rescaled through Eg. 1521 


change the amplitude of the cross-correlation function. This 
means that the cross-correlation function between two trac¬ 
ers does not have the bias 67^? and the correlation coeffi¬ 
cient = ^ab/(?a^b) is predicted to be smaller than unity. 
We therefore have two effects which can be used to constrain 
the relative velocity effect: (1) the BAO peak position and 
(2) the amplitude of the cross-correlation function relative 
to the auto-correlation functions. 


9.2 Constraining the relative velocity effect - 
BAO fits 

As mentioned in the last section, we have two effects which 
can be exploited to constrain the relative velocity effect, the 
shift in the BAO peak and the relative amplitudes of the 
individual correlation functions. While the shift in the BAO 
peak position can be considered as robust, there are effects 
other than the relative velocity bias which could change the 
amplitude. Any stochasticity, S 3 , in the galaxy density held, 
which is not correlated with the matter density, 5g — bg5m + 
Ss, and which does not correlate with the density held of 
the other survey, would lead to a reduction in the amplitude 
of the cross-correlation function. We will therefore perform 
multiple hts. First we will show the constraint on the relative 
velocity effect just using the BAO peak position, and then 
include the relative amplitudes of the correlation functions. 

We start with the BAO peak position. For this ht we 
convert the models of Eq.|3S]and[^into conhguration space 
and introduce additional polynomial terms to marginalize 
over the shape of the correlation functions, similar to the 
discussion in Section [S] 

?relvel(s) = B^Crelvel(os) -|- A{s) . (55) 

We also marginalize over the amplitude of the three correla- 
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Table 3. Summary of the fitting results for the relative velocity effect. The first two rows show the results when using only the BAO 
peak position to constrain bv before density field reconstruction. The last two rows present the same fit after reconstruction. The third 
and fourth rows list the fit parameters including the shape (and amplitudes) of the correlation functions. The fifth row provides the fit 
parameters to the correlation coefficient r^. All uncertainties are defined by the 68 % confidence levels. The fitting ranges are shown in 
the second column. 


fit condition 

fitting range 

by 

^CMASS 

WiggleZ 

^2 


pre-recon 



o.ot®i 

14 5+^-'^ 


BAO only 

50 - 180Mpc/fi 

-0.067 <bv < 0.010 

82.4/(78 - 16) 

BAO only 

50 - 180Mpc//i 

-0.31 <bv < 0.060 

1.0 

- 0.4 

84.0/(78 - 14) 

shape 

50 - 180Mpc//i 

-0.059 <bv < 0.096 

-2.01^® 

—0 7"^ ^ 

'-6.7 

89.4/(78 - 6) 

shape 

50 - 180 M.pc/h 

-0.12 <bv < 0.037 

1.0 

- 0.4 

94.7/(78 - 4) 


20 - 60Mpc//i 

—0.086 < bv < 0.062 

1.0 

- 0.4 

7.4/(9-4) 

post-recon 



'•^-8.3 

-OAtli 


BAO only 

50 - 180Mpc/fi 

-0.21 <bv < 0.02 

98.9/(78 - 16) 

BAO only 

50 - 180Mpc//i 

-0.22 <bv < 0.10 

1.0 

-0.4 

103.5/(78 - 14) 



Figure 18. The effect of the three relative velocity terms of Eg. 1481 
to the large scale correlation function. The black line shows the 
linear input correlation functions, while the magenta, red and 
blue lines include the bib^ and 6261 ) terms, respectively. The 
dashed lines represent the effect of the relative velocity terms 
(solid colored lines minus the black solid line). We fix bv = 0.1, 
bi = 2 and 62 = 1 - 


tion functions by giving each correlation function a separate 
bias parameter. Higher order terms for the cross-correlation 
function are always set by the bias of the auto-correlation 
functions. In total we have 14 free parameters {Bb, Be, Bw, 
Ab, Ac, Aw, ot, bv), where the polynomial terms A have 
three parameters each. Since we assume that the relative 
velocity is only present in CMASS, the relative velocity pa¬ 
rameter, bv, only affects the CMASS and cross-correlation 
function model. We perform fits where we additionally vary 
the parameter &2, but since this parameter is not well con¬ 
strained we often fix it to 1.0 for CMASS-BW and —0.4 for 
WiggleZ-BW (IYoo fc Seliakll201^ : iMarin et al.ll2013l ). The 
result of the fits are presented in the first two rows of Ta¬ 
ble [3] before reconstruction, and the last two rows after re¬ 
construction. Regardless of how we treat the parameter &2, 
we obtain constraints on which are consistent with zero. 



Figure 19. Distribution of the relative velocity bias bv obtained 
from the 480 mock catalogues using the “shape” and fitting 
methods (see Table [3] and text in Section l9.2l and l9At . 


This result is not surprising since the shift of the BAO peak 
due to the relative velocity effect is a shift to smaller scales. 
Our data, however, show a BAO peak at larger scales for 
CMASS-BW compared to WiggleZ-BW. 

9.3 Constraining the relative velocity effect - 
shape fits 

Next we fit the correlation functions without marginaliz¬ 
ing the relative amplitudes. In this case we include a bias 
parameter for CMASS-BW and WiggleZ-BW, but not for 
the cross-correlation function. The amplitude of the cross¬ 
correlation function is given by the product of the CMASS- 
BW and WiggleZ-BW bias parameters. Since this fit does 
not marginalize over the relative amplitudes and shape of 
the individual correlation functions, we also include the ve¬ 
locity dispersion parameter (t„ as 

Pg,Rna,i{k) = Pg{k)exp{-k^al/2), 

P^inAk) = Pg^{k) exp{-Aat/2). (56) 

Given the simplicity of our model, the parameter (j„ absorbs 
small-scale effects like non-linear structure formation which 
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are bias dependent. We therefore include three different cr„ 
parameters, one for each correlation function. The parame¬ 
ter (7„ introduces stochasticity on small scales which could 
mimic the relative velocity effect. We verify that this model 
yields a relative velocity bias consistent with zero when ap¬ 
plied to our mock catalogues. Figure [T^ shows the distribu¬ 
tion of obtained from the 480 mock catalogues. We only 
perform these hts pre-reconstruction, since we do not have a 
model for the post-reconstruction correlation function. The 
results are included in Table[3]with the label “shape”. While 
the constraints on become tighter compared to the “BAO 
only” fits, they are still consistent with zero. 

For the “shape” fit the correction factor of Eq. US] is 1.4 
and contributes significantly to our error budget. To avoid 
this additional source of error we can fit the correlation co¬ 
efficient instead of the individual correlation functions. 
This approach reduces the number of bins, which reduces 
the correction factor of Eq. 1151 The correlation coefficient 
should also be fairly independent of the underlying cosmo¬ 
logical model, since any effect common to the correlation 
functions cancels. However, the parameters bi, 62 and 
are degenerate when using . From Figure [TSl we can see 
that it might be possible to separately constrain 61 and 62 
given that 62 does behave very differently on small scales 
and around the BAO peak. However, while on small scales 
there is concern about the applicability of our model, on 
large scales the uncertainties are too large to exploit these 
effects. Thus we cannot vary all three parameters simultane¬ 
ously. We therefore fix the value of bi = 1.9 and 62 = 1.0 for 
CMASS-BW and bi = 1 and 62 = -0.4 for WiggleZ-BW. 
The term proportional to is usually significantly larger 
than the 61and 62&D terms, justifying to some extent our 
choice of fixing 61 and 62 (see Figure fTSll . 

The mock realizations are in good agreement with the 
expected value of the correlation coefficient above 20 Mpc/h. 
The scales above 60Mpc//i have large uncertainties and can 
be neglected for this fit, leading to the fitting range 20 - 
60 Mpc/h. We again verify that our model can reproduce 
a relative velocity bias of zero when applied to our mock 
catalogues. The distribution of maximum likelihood for 
the 480 mock realizations is included in Figure[T3(blue line). 
The best fitting parameters are included in Table |3| The 
relative velocity bias is again consistent with zero. 

So far the only constraint on the relative velocity bias, 

in the literature has been reported by IYoo fc SeliakI 
(l2013tl using the CMASS-DR9 power spectrum. They found 
bv < 0.033 at the 95% confidence level, consistent with our 
result. However, their constraint was found by fixing all cos¬ 
mological parameters to the Planck values, while our BAO- 
only constraint can be considered model independent. 


10 CONCLUSION 

We have investigated the galaxy clustering in the overlap 
region between the BOSS-CMASS and WiggleZ galaxy sur¬ 
veys. Having two galaxy samples in the same volume with 
different galaxy properties as well as survey selection effects 
presents a valuable opportunity to test for possible system¬ 
atic uncertainties in our analysis of the Baryon Acoustic 
Oscillation scale. We can summarize our results as follows: 

• We detect a BAO signal in both auto-correlation func¬ 


tions as well as the cross-correlation function of CMASS and 
WiggleZ using only the overlap region between the two sur¬ 
veys. The BAO detection in the cross-correlation function 
represents the first BAO detection in the cross-correlation 
function of two completely different galaxy surveys. After 
applying density field reconstruction we find distance con- 

-.fie! 

straints oi Dv^ = (1970 ± 45, 2132 ± 65, 2100 ± 200)Mpc 
for CMASS, the cross-correlation and WiggleZ, respectively. 
The three constraints are consistent with each other and 
with the distribution found in the mock realizations. The 
results are also robust against switching the displacement 
field of the two surveys during density field reconstruction. 
We therefore cannot see signs of systematic uncertainties. 

• We use our correlated mock realizations to determine 
the correlation between CMASS and WiggleZ. Using these 
correlations we derived a covariance matrix for the CMASS- 
DRll and our WiggleZ and cross-correlation constraints. 
While in the past the WiggleZ constraints have often been 
ignored when constraining cosmological models given the 
overlap (and hence correlation) with the CMASS results, 
our covariance matrix now allows one to make use of the 
WiggleZ information for cosmological constraints. Since the 
overlap region between the two surveys will not grow with fu¬ 
ture CMASS data releases, the covariance derived in this pa¬ 
per can easily be rescaled to obtain the covariance between 
our WiggleZ constraints and future CMASS data releases. 
We already provide a correlation matrix for the expected 
CMASS-DR12 results. 

• Using the measured correlation functions we test for the 
relative velocity effect, which is a possible source of system¬ 
atic uncertainty for BAO measurements. We perform various 
fits using the effect of the relative velocity bias on the BAO 
peak position as well as the relative amplitudes of the auto- 
and cross-correlation functions. We cannot detect any signs 
of a relative velocity bias. 
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